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ABSTRACT 

Using K-band imaging for 15 of the Canadian Network for Observational Cosmology (CNOCl) 
clusters we examine the near-infrared properties of moderate-redshift (0.19 < z < 0.55) galaxy clus- 
ters. We find that the number of K-band selected cluster galaxies within R500 (the Halo Occupation 
Number, HON) is well-correlated with the the cluster dynamical mass (M500) and X-ray Temperature 
(Tx)', however, the intrinsic scatter in these scaling relations is 37% and 46% respectively. Comparison 
with clusters in the local universe shows that the HON-M500 relation does not evolve significantly be- 
tween z = and z 0.3. This suggests that if dark matter halos are disrupted or undergo significant 
tidal-stripping in high-density regions as seen in numerical simulations, the stellar mass within the 
halos is tightly bound, and not removed during the process. The total K-band cluster light (L2oo.if) 
and K-band selected richness (parameterized by Bgcif) are also correlated with both the cluster T^; 
and M200- The total (intrinsic) scatter in the L2oo,i<"-M2oo and Bgc,_ff-M2oo relations are 43%(31%) 
and 35%(18%) respectively and indicates that for massive clusters both Ij2qo,k and ^gc,K can predict 
M200 with similar accuracy as T^, or optical richness (Bgc)- Examination of the mass-to- light 
ratios of the clusters shows that similar to local clusters, the K-band mass-to-light ratio is an in- 
creasing function of halo mass. Using the K-band mass-to-light ratios of the clusters, we apply the 
Oort technique and find flm,o = 0.22 ± 0.02, which agrees well with recent combined concordance 
cosmology parameters, but, similar to previous cluster studies, is on the low-density end of preferred 
values. 

Subject headings: cosmology: dark matter — large-scale structure of universe 
galaxies: clusters: photometry — fundamental parameters 



1. INTRODUCTION 

One of the principal objectives in modern cosmology is 
a good description of the physics that govern the forma- 
tion of structure from the galactic through to the super- 
cluster scale. Recent theoretical and observational work 
in this area has concentrated on understanding the rela- 
tionship between the dynamically dominant dark matter 
and the baryonic matter in the form of stars and gas. 
The major challenge has been that observables such as 
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the galaxy correlation function (e.g., Zehavi et al. 2005a, 
2005b; Eisenstein et al. 2005; Norberg et al. 2002, 2001), 
luminosity function (LF; e.g., Babbedge et al. 2006; II- 
bert et al. 2005; Dahlen et al. 2005; Wolf et al. 2003; 
Blanton et al. 2003; Nakamura et al. 2003; Kochanek 
et al. 2001; Beh et al. 2003; Cole et al. 2001), and the 
Lyman-alpha forest (e.g., McDonald et al. 2006; Kim et 
al. 2002; Croft et al. 2002) constrain the distribution 
of stellar mass and gas in the universe, yet dark matter 
is the dominant gravitational component, and numerical 
simulations are more effective at predicting its distribu- 
tion than that of baryonic mass.. 

Despite the difference between that which is easy to 
observe, and that which is easy to simulate, recent 
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combined N-body + Smoothed Particle Hydrodynamics 
(SPH) simulations have suggested that dissipation from 
the baryonic component is relatively unimportant in the 
process of galaxy formation because the force of grav- 
ity from galactic dark matter halos is overwhelmingly 
dominant (e.g., Weinberg et al. 2006; Nagai & Kravtsov 
2005). These results naturally explain why purely N- 
body simulations have been able to reproduce the corre- 
lation function (e.g., Colin et al. 1999; Berlind & Wein- 
berg 2002; Kravtsov et al. 2004; Tinker et al. 2005; 
Nagai & Kravtsov 2005; Tasitsiomi et al. 2004) and 
the number of satellite galaxies per halo, the Halo Oc- 
cupation Number (HON) or its probability distribution, 
the Halo Occupation Distribution, (HOD, e.g., Berlind & 
Weinberg 2002, see Cooray & Sheth 2002 for a review). 
Furthermore, full SPH simulations (e.g., Weinberg ct al. 
2006, 2004; Yoshikawa et al. 2001; Pearce ct al. 2001) 
and semi-analytic models (e.g., Berlind et al. 2005, 2003; 
Cole et al. 2000; Zhu et al. 2006) which can be directly 
compared to observations, have also begun to enjoy a 
great deal of success at reproducing the LF, correlation 
function, and HON. These recent improvements in the 
quality of simulations and semi-analytic models as well as 
new observational datasets such as the SDSS and 2dFRS 
represent a major breakthrough in our understanding of 
the relative distribution of baryonic and dark matter in 
the universe and consequently, the formation of large- 
scale structure. 

Detailed studies of galaxy clusters offer a complemen- 
tary approach to the correlation function, Lyman-alpha 
forest, and simulations for studying the relationship be- 
tween baryonic and dark matter. Clusters allow us to 
probe the distribution of baryons and dark matter in 
the most massive collapsed halos (which also happen to 
be the best-resolved objects in numerical simulations). 
The advantage of studying clusters is that the baryonic 
content in the form of stars and hot gas in the Intra- 
Cluster Medium (ICM) can be directly measured with 
observations, and in addition to this, cluster halos are 
sufficiently massive that their dark matter mass can be 
measured using either weak lensing. X-ray data, or the 
dynamics of cluster galaxies. This allows direct compar- 
ison between the baryonic and non-baryonic component 
on a halo-by-halo basis, whereas for galaxy-mass halos, 
the mass-to-light ratio (M/L), HON, or bias is usually 
measured statistically using either galaxy-galaxy lensing 
(e.g., Seljak ct al. 2005; Sheldon et al. 2004; Hoekstra et 
al. 2005, 2002), or by stacking galaxies and studying the 
dynamics of satellite galaxies (e.g., Conroy et al. 2005, 
Brainerd et al. 2003), although recently strong- lensing 
analysis of a sample of individual systems has been done 
(Treu et al. 2006). The advantage of being able to study 
individual halos is that not only can the correlation be- 
tween baryonic and dark matter be measured (via the 
HON or M/L ratio), but the relative scatter in the corre- 
lation can also be determined. Good constraints on the 
scatter are important because it is a direct measure of 
the stochasticity in the baryonic/dark matter bias. 

In practice, the best way to measure the baryonic con- 
tent of clusters in the form of stars is with detailed model- 
ing of the stellar populations using either high-resolution 
spectroscopy, or multi-band photometric observations. 
As these observations are often difficult to obtain for 
large samples of galaxies, K-band light has frequently 



been adopted as a cheap and efficient proxy for the stel- 
lar mass of galaxies. K-band light is neither strongly 
affected by dust, nor strongly enhanced by star bursts, 
and is therefore a good tracer of underlying stellar mass 
of a galaxy (e.g., Brinchmann & Ellis 2000; Gavazzi et al. 
1996; Rix & Rieke 1993). Furthermore, k-corrections are 
generally small, and fairly independent of galaxy type 
(Poggianti 1997; Mannucci et al. 2001). Therefore, a 
study of the relative abundances of K-band light and 
mass in clusters is a good probe of the relationship be- 
tween the cluster dark matter mass and the baryonic 
mass in the form of stars. 

In addition to being useful for determining the relation- 
ship between dark matter and baryonic matter, a good 
understanding of the slope and scatter of the cluster K- 
band luminosity-mass correlation and K-band richness- 
mass correlation will be extremely valuable in the era of 
high-yield cluster surveys such as the Red-sequence Clus- 
ter Survey 2 (RCS-2, Yee et al. 2007), the South Pole 
Telescope (SPT, Ruhl et al. 2004), and the Atacama 
Cosmology Telescope (ACT, Kosowsky 2006). These 
projects will use the evolution of the cluster mass func- 
tion, N(M, z) as a probe of cosmological parameters (see 
Mohr 2004 for a review); however, mass measurements 
for the thousands of clusters that will be detected in these 
surveys are unlikely to be available from traditional tech- 
niques such as Lj;, Tj;, weak- lensing, or dynamics and 
therefore new, observationally cheaper proxies for cluster 
mass will be required. Two candidates for this mass indi- 
cator are the integrated Sunyaev-Zeldovich Effect (SZE) 
?/-parameter (e.g., Motl et al. 2005), as well as the total 
cluster K-band luminosity or K-band selected richness 
(e.g., L04, Ramella et al. 2004, Rines et al. 2004). For 
the SZE surveys the y-parametcr has significant potential 
as a mass indicator (e.g., Hallman et al. 2006, Motl et 
al. 2005) as it is a measure of the mass of cluster baryons 
in the form of hot gas. For optical/IR surveys, K-band 
richness/luminosity is an attractive alternative because 
it is a measure of the total mass of baryons in the form 
of stars in cluster galaxies. Gladders et al. (2006) have 
already shown that optical richness works well as a clus- 
ter mass proxy, and using the optically-selected RCS-1 
survey and an empirical mass-richness calibration they 
determined cosmological parameters consistent with re- 
cent concordance values (e.g., Spergel et al. 2006), thus 
illustrating the potential for this technique. K-band light 
is a better tracer of stellar mass than optical light, and 
therefore it might be expected that the scatter in the K- 
band light/richness vs. mass relation will be smaller than 
the scatter in the optical light/richness vs. mass relation, 
or at least might have fewer catastrophic outliers. Using 
an observable with a smaller scatter, and fewer outliers 
could improve the accuracy in the measured cluster mass 
function (e.g., Lima & Hu 2005). 

Since the advent of the 2-Micron All Sky Survey 
(2MASS, Skrutskie et al. 2006) several studies of the 
relationship between K-band light and dark matter in 
local {z < 0.1) clusters have been done. Lin et al. (2003, 
2004; hereafter L04), Rines et al. (2004), Kochanek et 
al. (2003), and Ramella et al. (2004) aU show that 
the K-band selected number counts and total cluster K- 
band light are indeed correlated with the dark matter 
mass, and have rms scatters of ^ 40%. L04, Rines et 
al. (2004), and Ramella et al. (2004) also show that the 
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slope of the hn vs. M200 relation is shallower than unity 
for clusters, suggesting that the conversion of baryons 
into stellar mass proceeds less efficiently in higher-mass 
halos. Recently, Lin et al. (2006) presented an analysis 
of the K-band selected HON for a heterogeneous sample 
of 27 clusters at < 2 < 0.9 and found no significant 
evolution in the HON-Mass relation with rcdshift. In- 
terestingly, this result was different from their original 
analysis using a subsample of ^ 20 clusters at 0.2 < ^; < 
0.9 where they found the HON at fixed mass increases 
by a factor of ~ 2 at 2 > 0.2. 

The purpose of this paper is to extend the analysis of 
the local cluster K-band HON and M/L ratio to higher 
rcdshift using a well-defined sample of clusters with wide- 
field imaging and extensive spectroscopic data. Our sam- 
ple consists of 15 massive, X-ray selected clusters that 
were part of the Canadian Network for Observational 
Cosmology (CNOCl. Yce et al. 1996) project. The K- 
band imaging, LF and density profiles for the clusters 
were presented in the first paper in this series (Muzzin 
et al. 2006, hereafter Paper I). In addition to the K- 
band imaging, these clusters also have considerable opti- 
cal spectroscopy and g and r photometry both of which 
extend to R ^ R200 for each cluster. 

The structure of this paper is as follows: In §2 we 
briefly describe the dataset, and in §3 present updated 
values for the cluster dynamical masses and X-ray tem- 
peratures. In §4 we present the cluster HON and dis- 
cuss its redshift evolution. Section 5 shows that total K- 
band luminosity (L2oo,i<:) and K-band selected richness 
(parameterized by Bgc,/f ) are correlated with the cluster 
halo mass, and can potentially be used as cheap proxies 
for this quantity in cluster abundance surveys. In §6 we 
present the K-band M/L ratios for the CNOCl clusters 
and in §7 we present a measurement of flm using the 
Oort (1958) technique. We conclude with a summary in 
§8. When computing magnitudes and angular sizes we 
adopt an = 0.3, = 0.7, Hq = 70 km Mpc ~^ 
cosmology. 

2. DATA 

The CNOCl clusters are a set of 16 massive clusters 
with redshifts 0.17 < z < 0.54. Fifteen of the clusters 
were detected in the Einstein Medium Sensitivity Sur- 
vey (EMSS, Gioia et al. 1990). Abell 2390 was added as 
the 16th cluster and is also a massive cluster with strong 
X-ray emission. The clusters have extensive optical pho- 
tometry and spectroscopy which were obtained as part 
of the CNOCl project (Yce et al. 1996). Our sample 
is comprised of 15 of the 16 CNOCl clusters. The clus- 
ter MS0906-I-11 is omitted because it was shown to be a 
strong binary in redshift-space by Carlberg et al. (1996) 
and therefore the mass measurement for the cluster is 
unreliable. The optical photometric and spectroscopic 
data as well as the K-band photometric data used in this 
paper were already presented in Paper I. We refer to that 
paper as well as the CNOCl paper (Yee et al. 1996) for 
complete details of the observations, reductions and pho- 
tometry and below present only a quick overview of the 
data. 

2.1. Optical Photometry and Spectroscopy 

Gunn g and r band imaging data were obtained at 
the 3.6m Canada-Prance-Hawaii- Telescope (CFHT) us- 



ing the Multi-Object Spectrograph (MOS) camera. Pho- 
tometry was performed on these data using the Picture 
Processing Package (PPP, Yee 1991). The photometry 
reaches a 5a depth of ^ M* -|- 3 in both the g and r 
bands. The CNOCl collaboration also obtained > 2500 
spectroscopic redshifts in the fields of the 15 clusters us- 
ing the CFHTMOS. Of these, approximately one-half are 
cluster members. The spectroscopy is sparsely sampled, 
but complete to a depth of K* -I- 2 for all but the two 
highest-redshift clusters (MS0016-hl6 and MS0451-03). 
The spectroscopy for those clusters is complete to ~ K* 
+ 1. 

2.2. Near Infrared Photometry 

K-band imaging for 14 of the 15 clusters was obtained 
at the Kitt Peak National Observatory (KPNO) 2.1m 
telescope using the Ohio State / NOAO Infrared Imag- 
ing Spectrograph (ONIS). Observations were taken in a 
Mauna Kea filter set version of the K^ filter (Tokunaga 
K), which is nearly identical to the 2MASS Kg filter. For 
our analysis we treat the Tokunaga K filter as the Kg- 
band and hereafter refer to it as the "K-band" . K^-band 
imaging of MS0440+02 was obtained using the PISCES 
camera on the Steward Observatory 90" telescope. Pho- 
tometry for the K-band data was also performed using 
PPP. The K-band imaging is complete to a 5<j depth of 
~ K* -I- 2 for all clusters, except the two highest-redshift 
clusters where it is complete to ~ K* -t- 1 (see Table 1 of 
Paper I for a summary of observations). 

3. CLUSTER PHYSICAL PARAMETERS 

The masses of the CNOCl clusters are without a doubt 
the most well-studied for clusters at moderate rcdshift. 
There are numerous measurements using X-ray temper- 
atures {Tx, Hicks et al. 2006; Lewis et al. 1999; Henry 
2000; Mushotzky & Scharf 1997), the dynamics of cluster 
galaxies (van der Marel et al. 2000; Borgani et al. 1999; 
Carlberg et al. 1996; 1997; Diaferio et al. 2005), and 
both strong (Fahlman et al. 1994; Luppino & Gioia 1995; 
Pierre et al. 1996, Wu 2000) and weak lensing (Allen 
1998; Hoekstra et al. 1998; Small et al 1995; 1997). In 
this section we discuss the masses and dynamical radii 
we adopt for comparing with the K-band properties. 

3.1. Dynamical Masses 

The first study of the dynamics of the CNOCl clus- 
ters was done by Carlberg et al. (1996) who measured 
the line-of-sight velocity dispersion (ci) for each clus- 
ter. These velocity dispersions were used to calculate the 
virial mass (M^j^), the radius at which the cluster mass 
density exceeds the critical density of the universe by a 
factor of 200 (R200), and the mass contained within that 
radius (M200), by assuming the cluster had an isotropic 
velocity ellipsoid. In subsequent work, Carlberg et al. 
(1997) and van der Marel et al. (2000) examined in detail 
the velocity dispersion profiles, and velocity anisotropy 
of the clusters, van der Marel et al. (2000) verified that 
they are compatible with an isotropic velocity ellipsoid 
and that the cluster mass profile was consistent with an 
NFW (Navarro et al. 1997) profile. The cluster velocity 
dispersions were also determined by Borgani et al. (1999) 
using different background subtraction techniques. They 
found that the different techniques provide velocity dis- 
persions that are self-similar to 10%, and that the 
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velocity dispersions they derive arc also consistent to ~ 
10% with the Carlbcrg et al. (1997) velocity dispersions. 

Unfortunately, the M200 and R200 values determined 
from these studies are now out of date because they were 
computed assuming cosmological parameters which are 
different from recent concordance values (e.g., Spergel et 
al. 2006). Fortunately, the velocity dispersions do not 
depend on cosmology and, given that the clusters are 
consistent with an isotropic velocity ellipsoid, we have 
recomputed both R200 and M200 with a Q,n = 0.3, f^A = 
0.7, Ho = 70 km s~-^ Mpc~-^ cosmology using the equa- 
tions, 

■^200 = 7777773:: (1) 



and. 



M200 = 37ri?2oo • 200pc 



(2) 



where ii{z) is the Hubble constant at redshift z, and Pc 
is the critical density for a flat universe. Equation 1 can 
be derived by assuming M200 ~ Myir and equating the 
virial theorem (M„jr = ^(^fRvir) and Equation 2. For 
all clusters the ci values determined by Carlbcrg et al. 
(1997) are used to compute M200 and R2oo- Errors in 
M200 are calculated by propagating the errors in ai in 
the standard way. The new values of R200 and M200 are 
on average ^ 1.6 times larger than the Carlberg et al. 
(1997) values. We list the updated values in Table 1. 

3.2. X-ray Temperatures 

X-ray temperatures have been measured for the 
CNOCl clusters by several authors. The clusters wc 
originally discovered in the Einstein Medium Sensitivity 
Survey (EMSS, Gioia et al. 1990); however, only X-ray 
fluxes were computed using the original data. Subse- 
quently, Lewis et al. (1999) measured T^, for 13 of the 15 
clusters in our sample using data from the ROSAT satel- 
lite. They found that the average cluster masses deter- 
mined from Tj; were in good agreement with the updated 
dynamical masses, even though individual clusters could 
have discrepancies as large factor of 2. Recently, 
Hicks et al. (2006) used archival Chandra Advanced 
CCD Imaging Spectrometer (ACIS) data to determine 
accurate T^, values for 13 of the 15 clusters in our sample. 
They use single and double-/3-model fits to determine 
the temperatures. Several clusters (A2390, MS0839-F29, 
MS1358+62, MS1455-h22) are also corrected for signif- 
icant cooling flows. They find that the X-ray masses 
they compute agree well with the dynamical masses de- 
termined with the velocity dispersions of Carlberg et 
al. (1997) and Borgani et al. (1999). These new T^, 
values are superior to the older data and for compar- 
isons with the HON (§4), L2oo,K (§5.1), and Bgc,K (§5.2) 
we use these values. There are two clusters without a 
Chandra observation (MS1224-I-20, MS1231-M5). For 
MS1224-I-20 we use the temperature listed in Lewis et al. 
(1999). MS1231-I-15 has no X-ray temperature available 
so we assume a temperature of 6 keV, which is appropri- 
ate given the X-ray luminosity (Yee & EUingson 2003). 
We list the X-ray temperatures and their associated er- 
rors in Table 1. 

4. THE HALO OCCUPATION NUMBER 



Clusters are massive dark matter halos, and therefore 
computing their HON simply requires counting the total 
number of cluster members within some dynamical ra- 
dius, typically R2oo- The extensive spectroscopy for the 
CNOCl clusters could be used to differentiate between 
field and cluster galaxies; however, we prefer to mea- 
sure the HON using statistical background subtraction. 
Using spectroscopic data to separate field and cluster 
galaxies is preferable when computing cluster properties 
that depend on the luminosity of cluster members (e.g., 
LFs, total luminosity §5.1, or M/L ratios §6.1) because 
these calculations require information on the member- 
ship of individual galaxies to compute distance moduli 
and k-corrections. However, unless the spectroscopy of 
a cluster field is quite complete, statistical background 
subtraction is better-suited for counting the overdensity 
of cluster galaxies. Furthermore, using statistical back- 
ground subtraction provides a consistent technique for 
each cluster and avoids any biases from cluster-to-cluster 
that might be caused by poor determination of the spec- 
troscopic selection function. 

The background counts are measured from our own 
K-band imaging survey of CN0C2 fields (Yee et al. 
2000). The data was taken using the PISCES camera 
on the Steward Observatory 90" telescope and are re- 
duced and photometered using the same techniques as 
for the CNOCl dataset. The imaging consists of 23 (8' 
X 8') fields and covers a total area of 0.26 square degrees. 
Complete details of the observations, data reduction, and 
photometry will be presented in a future paper (Lin et 
al. in preparation). Figure 1 shows background counts 
from the CN0C2 fields with the number counts from the 
K20 Survey (Cimatti et al. 2002) and the MUNICS Sur- 
vey (Drory et al. 2001) overplotted for comparison. The 
CN0C2 galaxy counts agree well with the counts from 
these surveys as well as various other K-band surveys in 
the literature (e.g., Elston et al. 2006; Maihara et al. 
2001; Saracco et al. 2001; Gardner et al. 1993). 

We measure the HON for the clusters within the ra- 
dius R500 (the radius at which the cluster mass density 
exceeds the critical density of the universe by a factor 
of 500), rather than the more typical R2oo- It would 
be preferable to use R200, as it is similar to the clus- 
ter virial radius; however, some of the clusters are very 
sparsely sampled at R200, and therefore require large cor- 
rections to compensate for the poor sampling. These 
corrections propagate as large uncertainties in the total 
number of galaxies. Conversely, the coverage at R500 is 
complete for all clusters (except MS1455+22) and there- 
fore those measurements have much smaller uncertainties 
and should also be more robust. The cluster R500 is es- 
timated from the R200 values assuming an NFW profile 
with c = 5. Hereafter we refer to the HON within R500 
as N500 and the HON within R200 as N2oo- 

Thus far, the largest study of the cluster HON is the 
work of L04 who measured N500 and N200 for 93, X-ray 
selected, ^; < 0.1 clusters using 2MASS data. In order to 
make a direct comparison between our sample and their 
sample, wc measure N500 using the same limiting mag- 
nitude as L04. Their technique for measuring the HON 
involves determining the Schechter parameters K* and 
0* for each cluster and then integrating the Schechter 
function to an absolute magnitude limit of M^ = -21 to 
find the total cluster counts. For their ensemble of 93 
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clusters, they measured K* — -24.02 ± 0.02. Therefore, 
statistically speaking, their N500 values are measured at 
K* + 3, although the depth varies from cluster-to-cluster. 
Our data does not reach K* -|- 3 for most clusters, there- 
fore galaxies are counted to K* -|- 1 and then these counts 
are extrapolated to K* -I- 3. 

K* evolves with redshift; however, we showed in Paper 
I that the evolution tightly follows the luminosity evolu- 
tion of a single-burst population formed at high-rcdshift. 
The single-burst model from Paper I is used to infer K* at 
the redshift of each cluster. The background subtracted, 
K < K* + 1 counts are scaled to the appropriate counts 
for K* + 3 by extrapolating a Schechter function with a 
faint-end slope of a = -0.9. This value of a is consistent 
with the faint-end slope measured for the CNOCl cluster 
LF in Paper I (a = -0.84 ± 0.08) and is similar to the 
faint-end slope measured for the L04 data (a = -0.84 ± 
0.02). 

We plot N500 vs. M500 for the clusters in the left panel 
of Figure 2, and N500 vs. T^: in right panel of Figure 2. 
In both cases there is a good correlation between the two 
parameters. The Spearman rank-correlation coefficients 
are 0.79 and 0.77 respectively, implying probabilities of 
7.1 X 10~^, 8.5 X 10"'' that the data are uncorrelated. 

One cluster which has an extremely small N500 for its 
mass is MS1455-I-22. This cluster is also a significant 
outlier in the M2oo-Bgc relation (Bgc is a measure of 
cluster richness, see §5.2) from Yee & EUingson (2003, 
hereafter YE03). They suggest that the most logical in- 
terpretation is either that the mass has been significantly 
overestimated or else that the cluster is in a very ad- 
vanced state of evolution, different from most clusters in 
this redshift range. Therefore, MS1455-f22 is excluded 
when we fit both M500 vs. N500 and and T^; vs. N500 
and is plotted in Figure 2 as an open triangle. Using a 
X^-minimization technique for errors in both parameters 
(Press ct al. 1992) we find the best-fit relation for the 
N500 - M500 correlation is 

Log{M5oo) = (1.40 ± 0.22)Lo5(iV50o) + (11-58 ± 0.54), 

and has a reduced-x of 1.60. For the N500 - T^, the 
best-fit relation is 

Log{T^) = (0.77±0.11)Lo5(iV50o) + (-1.01±0.27), (4) 

and has a reduced-x^ of 2.46. The correlations have 

Root-Mean Squared (rms) scatters of 51% and 50% re- 
spectively. For the remainder of the paper, the percent- 
age rms scatter in all correlations is computed using the 
same method as YE03. The rms deviation from the fit 
in terms of the dependent variable (in this case, M500) is 
determined and then compared to its mean value from 
the entire sample. 

The scatter in these correlations is larger than the 35% 
scatter in the local scaling relations measured by L04. 
However, the random errors in our data are larger than 
theirs, and clearly a portion of the scatter comes from 
these measurement errors. Given that the total scat- 
ter is the result of both the measurement errors in M500 
and N500 as well as some intrinsic scatter, we estimate 
the intrinsic scatter by assuming that the measurement 
scatter is independent of the intrinsic scatter and then 
subtracting the mean measurement scatter of the depen- 
dent variable in quadrature from the total scatter. Using 



this simplistic method we estimate that the intrinsic scat- 
ter in N500 - M500 is 37% and for N500 - T^, it is 46%. 

We compare the N500 values from the CNOCl clus- 
ters with the L04 values in the top panel of Figure 3. 
L04 measured N500 by integrating the cluster LFs over 
an NFW spatial distribution and therefore it is the N500 
within a spherical volume. Our N500 is measured us- 
ing statistical background subtraction and is therefore 
the number of galaxies in a cylinder of radius Rsoo- We 
convert our N500 in cylinders to N500 in spheres by multi- 
plying by a deprojection constant of 0.791, the value for 
the percentage difference in enclosed number of galaxies 
between spheres and cylinders at R500 for an NFW pro- 
file with c = 5. The dash-dot blue line in Figure 3 is 
the best-fit relation for the L04 clusters, and the solid 
red line is the best-fit relation for the CNOCl clusters. 
The slope of the CNOCl relations (M500 cx N|?oo^=^"-") 
is slightly shallower than the L04 relation that includes 
the BCGs (M200 oc N°p^(f=^° °''), but is consistent within 
the errors. It is also consistent with the scaling relation 
found by Rines et al. (2004) for the CAIRNS clusters 
(M200 oc ^200^^^^)'^ however, their relation is measured 
using N200 and M200 and a different luminosity cut. The 
overall normalization of the CNOCl N500 is nearly identi- 
cal to the L04 normalization. As a comparison between 
the samples we plot the measured N500 values divided 
by those predicted from the local M500-N500 relation in 
the bottom panel of Figure 3 . Figure 3 shows that 
the HON - mass correlation measured by L04 does not 
evolve between z = and z ^ 0.3. The higher redshift 
CNOCl clusters have N500 values which are only 4% ± 
11% smaller than would be predicted by the local rela- 
tion. The mean(median) N 500/^ 500 jitiocai is 0.96(0.90) 
± 0.11, and is consistent with no evolution in the HON 
with redshift. As a comparison, the L04 clusters have a 
mean(median) ^500 f^^boa, fttiocai of 1.04(1.01) ± 0.04. 

The non-growth of the HON at a fixed mass with in- 
creasing redshift is different from the results of L04. In 
addition to their local clusters, L04 also computed the 
HON for a subset of clusters with 0.1 < ^; < 0.9 (K-band 
imaging from De Propris et al. 1999) which had avail- 
able X-ray temperatures and compared it to the local 
clusters. Their comparison showed that the HON was 
roughly a factor of 2 larger in the z > 0.1 clusters. Prom 
this they suggested that half of all cluster galaxies may 
be removed (via either merging or disruption) between 
z ^ 0.6 (all but two of their clusters are z < 0.6) and 
z = in order to match the local HON. 

It is unclear why the conclusions of our studies are so 
different. We did note in Paper I that it is challenging 
to reconcile the idea of a strong evolution of the cluster 
HON between z ~ 0.6 and z = and the passive evo- 
lution of the K-band cluster LF. Paper I showed that 
the evolution of the cluster K* tightly follows the predic- 
tion of a passive evolution model between z ~ 0.5 and 
z = and those results reaffirmed the findings of the 
pioneering cluster K-band luminosity function study of 
de Propris et al. (1999), as well as recent K-band LFs of 
high-redshift clusters (e.g.. Toft et al. 2003; Kodama & 
Bower 2003; Ellis & Jones 2002; Strazzullo et al. 2006). 
Interestingly, all of the 2; > 0.1 clusters used by L04 are 
part of the de Propris et al. (1999) sample. If 50% of 
cluster galaxies are either merged, or disrupted by tidal 
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forces between z ^ 0.6 and z = 0, then it seems unlikely 
that the evolution of K* could follow a passive evolution 
model, unless the probability of being merged or tidally 
disrupted was independent of galaxy K-band luminosity 
(i.e., stellar mass). While mass-independence may be 
plausible in the case of mergers, more massive galaxies 
are more difficult to disrupt. If a significant number of 
disruptions do occur, they should preferentially destroy 
galaxies on the faint-end of the luminosity function and 
therefore will change its overall shape. 

L04 appeal to the N-body simulations of Kravtsov et 
al. (2004) for an explanation of the strong redshift evo- 
lution of the HON. Kravtsov ct al. (2004) show that for 
dark matter halos, the HON at a given mass becomes 
approximately a factor of three larger from .2 = to ^ = 
5. They show that the dark matter HON at a given 
mass is smaller at lower redshift because low mass ha- 
los merge with high mass halos as well as become tidally 
truncated or disrupted by massive halos in high-density 
regions such as clusters. 

More recent simulations which incorporate both dark 
matter as well as baryonic matter using Smoothed Parti- 
cle Hydrodynamics (SPH, e.g., Nagai & Kravtsov 2005, 
Weinberg et al. 2006) show that when the baryonic com- 
ponent is included, the interpretation of the HON de- 
rived by counting galaxies becomes more subtle. Nagai 
& Kravtsov (2005) simulated a set of 8 clusters and found 
that as subhalos enter the cluster R200 they are stripped 
of ^30% of their total dark matter mass. Once they fall 
into the cluster they are further stripped until ~70% of 
their original dark matter mass is removed. Their result 
confirms the earlier N-body simulations which showed 
the dark matter HON is depleted with decreasing red- 
shift, especially in high-density environments; however, 
the SPH simulations also show that primarily the dark 
matter is stripped from the halos and very little of the 
baryonic mass (in the form of stars and gas) is lost dur- 
ing this process. The baryonic matter remains because 
it tends to cool, and fall to the center of the subhalo and 
is amongst the most tightly bound mass in the subhalo. 

These simulations suggest that the HON, when se- 
lected by counting the number of massive dark matter 
subhalos is quickly depleted in high density regions such 
as clusters because of tidal stripping and merging of 
the subhalos. However, when the HON is selected by 
the baryonic mass of halos (i.e., using K-band number 
counts) there should be little evolution in the HON at a 
given halo mass, because the baryonic component of sub- 
halos is tightly bound, and will only be depleted in the 
case of full disruption. This interpretation is consistent 
with our data and is also consistent with the observa- 
tion that both K* (Paper I, L04, de Propris et al. 1999) 
and Q (L04) for cluster galaxies are basically identical 
across more than an order-of magnitude in cluster mass 
(Paper I, L04, de Propris et al. 1999). Given that the ef- 
ficiency of tidal stripping is correlated with cluster mass, 
it is hard to understand how the K-band LF of cluster 
galaxies could be independent of cluster mass unless the 
baryonic component is unaffected in the process. 

In summary, we find that the K-band selected HON 
of clusters at a given mass does not evolve significantly 
between 2; = and z ^ 0.3. This result is consistent with 
the more recent work of Lin et al. (2006) who measured 
the HON within R2000, Riooo and R700 for a sample of 



27 clusters < z < 0.9 (a few of which are included 
in L04) and also find no-evolution in the HON at fixed 
mass. We argue that the non-evolution of the HON is 
naturally explained by recent SPH simulations, which 
show that while dark matter halos encounter a signif- 
icant amount of tidal stripping in high-density regions 
such as clusters, the baryonic component within remains 
mostly intact. Furthermore, this interpretation recon- 
ciles the strong tidal stripping with the purely passive 
evolution of the cluster K-band LF, and its invariance 
across the cluster mass spectrum. 

Our data cannot rule out the possibility that the simu- 
lations may over-exaggerate the tidal stripping and that 
similar to the baryonic component, the dark matter 
mass-selected HON at a fixed cluster mass also does not 
evolve within the same redshift range. A useful way to 
test this would be to compare the M/L ratios of clus- 
ter galaxies and field galaxies at large radii and look for 
evidence of truncated halos within the cluster popula- 
tion. Unfortunately, this is a challenging task, because 
there are few tracers of the dark matter potential at large 
radii for galaxies. There is now evidence from galaxy- 
galaxy lensing studies of massive clusters that supports 
the tidal-stripping scenario (Halkola et al. 2006, Natara- 
jan et al. 2002); however, in lower-mass clusters (M ~ 
10^^-10^^ Mq), the same trend is not observed (Man- 
dclbaum et al. 2006). A more comprehensive study of 
galaxy-galaxy lensing in clusters in the range of M ~ 
10^^-10^'' Mq would be particularly useful for testing this 
interpretation. 

5. K-BAND LIGHT AND RICHNESS AS AN INDICATOR OF 
CLUSTER MASS 

In this section we measure the total K-band luminos- 
ity within R200 for the clusters (L2oo,a:) and examine its 
correlation with cluster mass. We also determine the 
correlation between K-band selected richness and cluster 
mass. 

5.1. Tot,al K-Band Luminosity 

J-'200,K for a galaxy cluster is defined as the sum of the 
K-band luminosity of all cluster galaxies within R200 to 
a fixed absolute magnitude. In principle the measure- 
ment is straightforward; however, in practice, data are 
never homogeneous and additional corrections must be 
applied. In the case of the CNOCl clusters, the clus- 
ters lie at a fairly large range of redshifts and therefore 
the light from each galaxy must be both k-corrected and 
evolution corrected to a common redshift. The need for 
these corrections means that the membership and ap- 
proximate spectral-type for each galaxy in the cluster 
field must be determined. This is not as precise us- 
ing statistical background subtraction and therefore we 
make use of the extensive spectroscopic catalogues. The 
spectroscopy is complete to K* -|- 1 for all clusters, but 
is sparsely sampled. Therefore, a spectroscopic selection 
function is used to correct for the sampling. The method 
for determining the spectroscopic selection function was 
developed by the CNOCl (Yee et al. 1996) & CN0C2 
(Yee et al. 2000) collaborations and is discussed in de- 
tail in those papers as well as Paper I. Similar to the 
analysis in Paper I, we ignore the "secondary" selection 
effects (color, position, and z) and use only the mag- 
nitude weights, which are the overwhelmingly dominant 



Muzzin et al. 



7 



10^ 



10'' 



in 
p 



1 102 



lOi 



10" 



-n i I i i i p 

• CNOC2IR 

■ MUNICS - Drory et al. (2001) 

▲ K20 - Cimatti et al. (2002) 



A A ik A 



I I 




A. "A 



12 



14 



16 



18 



20 



22 



K(Mag) 



Fig. 1. — Red Dots: Number of galaxies per square degree per 0.5 magnitudes as a function of magnitude determined from 0.26 square 
degrees of CNOC2 K-band imaging data. Green Triangles: Number counts from the K20 Survey. Blue Squares: Number counts from the 
MUNICS Survey. 
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Fig. 2. — Left Panel: Plot of M500 vs. N500 for the CNOCl clusters. The solid line is the best-fit linear relation and has an rms scatter 
of 50%. Right Panel: Plot of T^, vs. N500 for the same clusters. The solid line is the best-fit linear relation and has an rms scatter of 51%. 
MS1455-I-22 is plotted as an open triangle. 

directly compared to local studies. 

The Ij2oo.k for some clusters must also be corrected for 
the incomplete coverage of R2oo- We noted that for the 
HON (§4) the coverage corrections propagate to large er- 
rors in the HON when it is computed within R200 and 
because of that, the HON was measured within Rsoo- 
The incomplete radial coverage is less problematic when 
computing Ij2oo,k because much of the total cluster lu- 
minosity comes from the BCG (~ 5 - 30%). Despite the 
fact that the number of cluster galaxies at R500 < R < 
R200 is large, the most luminous galaxies tend to reside 
in the cluster core, and therefore galaxies at R500 < R < 
R200 contribute less to the total luminosity than to the 
HON. 



selection bias (Yee et al. 1996). 

The k-correction for each galaxy is taken from the 
models of Poggianti (1997). The k-corrections depend 
only mildly on spectral- type and the color/spectral- type 
model discussed in Paper I is used to estimate a spectral- 
type for each galaxy. Also, the LPs from Paper I showed 
that the luminosity evolution in the clusters is passive, 
but significant (0.35 ± 0.06 magnitudes from z = to z 
^ 0.5). Therefore, an evolution correction from the Pog- 
gianti (1997) models is applied to each galaxy, where we 
have transformed the Poggianti (1997) evolution correc- 
tions to our cosmology (see Paper I, §6.1.2). By applying 
both k and evolutionary corrections the measured L2oo,if 
is corrected to z = 0, and therefore these values can be 
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Fig. 3. — Top Panel: Plot of N500 vs. M500 for the CNOCl clusters (large red points) and the L04 clusters (small blue points). 
MS1455-I-22 is plotted as an open triangle. The solid red line is the best-fit relation for the CNOCl clusters, and the dash-dot blue line 
is the best-fit relation for the L04 clusters. The slope is the CNOCl correlation is slightly shallower, but consistent with the L04 slope. 
Bottom Panel: Plot of the measured N500 divided by the value predicted by the L04 scaling relation as a function of redshift for the L04 
clusters (blue points) and CNOCl clusters (red points). The CNOCl normalization is consistent with the L04 normalization suggesting no 
evolution in the HON with redshift. 



The clusters without complete coverage have been ob- 
served in a strip through the cluster core; therefore, the 
coverage problems are corrected by computing the lumi- 
nosity in circular shells. The total luminosity of a shell 
is multiplied by the ratio of the total area of the shell, 
to the total area of the shell with imaging/spectroscopic 
data. Several of the clusters {MS0016-hl6, MS0451-02, 
MS1006-hl2, MS1008-12, and MS1455-h22) have no cov- 
erage for a few of the outermost shells. For these clusters 
a correction based on the profile of clusters which have 
the best coverage in the outer regions is applied. This 
correction is small in four of the five cases (~ 5-10%), 
again because the majority of the cluster light comes 



from the core. MS1455-I-22 has a very large R200 (which 
is likely to be overestimated, see §4) and therefore the 
correction is much larger (~ 60%). 

Once the selection function, k-correction and evolution 
corrections have been applied, and the sum of the lumi- 
nosity of galaxies brighter than K* -I- 1 has been tabu- 
lated, the "total" luminosity of the cluster is finalized by 
extrapolating the LF to a fixed absolute magnitude. We 
choose to extrapolate to Mk = -21 (i.e., K* -I- 3 at z 
= 0), the same limiting absolute magnitude used by L04 
and and effectively the same limiting magnitude of the 
HON. We evolution-correct the combined LF from all 15 
clusters to z = (K* = -24.14 ± 0.15 and a = -0.84 ± 
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0.08, sec Paper I) and use this for extrapolation. This 
LF is consistent with the average LF found by L04 (K* 
= -24.02 ± 0.02 and a = -0.84 ± 0.02) and because our 
clusters have been evolution corrected to 2; = 0, extrap- 
olating the luminosity to = -21 results in a L2oo,if 
which can be directly compared between the samples. 
Before the correction is applied, the contribution from 
the BCGs is removed from the total cluster light. The 
BCGs contribute a significant, but variable, fraction of 
the total cluster Ught (~ 5-30% at K* + 1) and do not 
obey a Schechter function. Including them as part of 
the extrapolation of the Schechter function would over- 
estimate the total cluster light. The light from the BCG 
galax;y is re-added to the total once the extrapolation has 
been done. Lastly, a bootstrap error for each cluster is 
computed by performing the entire analysis using 300 re- 
samplings with replacement from the data. Table 2 lists 
the final lj2m,K values as well as the bootstrap errors. 

We examine the correlation between the cluster mass 
and light by plotting Log(M2oo) vs. Log(L2oo,if ) in the 
first panel of Figure 4. The best-fit linear relation is 

Log{M2oo) = {1.20±0.16)Log{L2oo,K) - (0.95 ± 2.21), 

(5) 

and the fit has a reduced-^ of 2.65, where again we have 
ignored MS1455+22 because it is likely to have an incor- 
rect M200- The rms scatter in the M200 - ^200,k relation 
is 43% and from this we estimate that the intrinsic scat- 
ter is 31%. The scatter in the correlation is about a factor 
of 1.3 larger than the 34% total, and 24% intrinsic scat- 
ter measured by L04 for local clusters. If we invert the 
axes, we find L2oo,x « M^qq'^"'^"'^^. This agrees with the 
scaling relation for local clusters with the BCG included 
which is L200.A: oc M^qq^'^'''''^ (L04), and somewhat less 
well with the scaling relation for local groups which is 
UooM « M°ao^"-°^ (Ramella et al. 2004). This shows 
that the slope, while slightly shallower for groups, does 
not change significantly over 2 orders of magnitude in 
cluster mass. Furthermore, the reasonable agreement be- 
tween the slope and intrinsic scatter in the correlations 
at different redshifts suggests that whatever physical pro- 
cesses are responsible for building in the correlation have 
largely taken place by z ~ 0.3, and that there is little 
change in the L200 - M200 relation from z ^ 0.3 to 2; = 
0, besides possibly a small decrease in the scatter. 

In the second panel of Figure 4 we plot the correlation 
between the cluster X-ray temperature and L2oo,if- The 
best fit relation with MS1455-I-22 excluded is 

Log{T^) = (0.77 ± Om)Log{L2oo,K) + (-9.45 ± 1.07), 

(6) 

and has a reduced-x of 2.35. The rms scatter in the re- 
lation is 24% and this corresponds to an intrinsic scatter 
of 14%, notably better than the M200 - ^200, k relation. 
As the CNOCl clusters are all massive, relaxed X-ray 
systems, this may not be surprising. Using r-band data, 
YE03 found the correlation between Bgc and T^; for the 
CNOCl clusters to have a smaller scatter (21%) than 
the Bgc - M200 relation (31%), although they used the 
Tx values determined by Lewis et al. (1999) rather than 
the Hicks et al. (2006) values. 

If we compare the scatter in the M200 vs. L2oo,k rela- 
tion to the scatter in the correlations of M200 vs. T^;, L^,, 
and Bgc (optical) determined by YE03 for these clusters, 
we find that the total K-band light is approximately as 



accurate as those parameters in inferring the dynamical 
mass. YE03 showed the scatter in the M200 vs. T^;, 
La;, and Bgc (optical) for these clusters was 29%, 35%, 
and 31% respectively. This demonstrates that L2oo,a" is 
as good as, but no more accurate than traditional mass 
indicators at inferring the cluster dynamical mass. It ap- 
pears the fact that K-band light is a fairly clean tracer 
of an individual galaxy's stellar mass does not result in 
a smaller scatter for the L2oo,k - M200 scaling relation, 
suggesting that the amount of scatter measured in opti- 
cal scaling relations is a real scatter in relative amounts 
of stellar mass and dark matter mass in clusters, and not 
enhanced due to a variance in star-formation properties 
or stellar populations from cluster-to-cluster. 

It is possible that the scatter we measure is enhanced 
by systematics in the data analysis caused by the need 
for k-corrections, evolution corrections, spectroscopic se- 
lection functions, and coverage corrections; however, the 
level of scatter we find is of order that found by L04 in lo- 
cal clusters, which do not suffer from those uncertainties. 
Given that there is no reason to expect the scatter in the 
scaling relations to become smaller at higher-redshift, 
it would suggest that these corrections have been ap- 
plied properly and do not significantly increase the over- 
all scatter. 

The correlations from Figure 4 are pleasing in the sense 
that total K-band light is observationally much cheaper 
to measure than L^; , T^; or M200 (via a velocity dispersion 
or weak lensing) for a given cluster. For clusters in this 
redshift range, a wide-field NIR camera on a 4m-class 
telescope can achieve a limiting magnitude of < K* 
-|- 1 with an integration time of only a few minutes. This 
means followup of many candidates in the generation of 
high-yield cluster cosmology projects is easily achievable. 
Unfortunately, a major concern for cosmology projects is 
that although L2oo,k appears to be a good proxy for M200 
or Tj; it cannot be used in practice because it requires an 
a priori knowledge of R2oo- In the next section we show 
that the NIR-selected richness in a fixed aperture is also 
well-correlated with the cluster mass and can be used as 
a mass proxy. 



5.2. Bgc as an Indicator of Cluster Mass 

The good correlation between L2ao,K and M200 found 
for the CNOCl clusters, and by other authors (e.g. Lin 
et al. 2003; L04; Rines et al. 2004; Ramella et al. 2004) 
demonstrates that K-band light can be used as an ef- 
ficient estimator of the cluster halo mass. Of course, 
measuring L2oo,if requires knowledge of R20O) which is 
generally an unknown, whereas the cluster richness, pa- 
rameterized by Bgc does not require R2oo- Bgc is the 
amplitude of the 3-dimensional, spatial correlation func- 
tion between the cluster center and the cluster galax- 
ies. If the shape of the 2-dimensional angular correlation 
function is assumed (i.e., w{9) oc d'--^ 'K^ with 7 ~ 1.8), 
then its amplitude can be measured by counting galax- 
ies in a fixed aperture around the cluster center. Bgc is 
then measured by deprojecting the angular correlation 
function via (^(r) oc r~'^) and scaling it by a luminosity 
function. A full derivation and motivation for Bgc is pre- 
sented in Longair & Seldner (1979). The Bgc parameter 
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Fig. 4. — Left Panel: Plot of Log(M2oo) vs. Log(L2oo,if ) for the CNOCl clusters. The solid line is the best-fit linear relation. The rms 
scatter around the relation is 43% (31% intrinsic). Right Panel: Plot of T^, vs. Log(L2oo,/f )• The solid line is the best-fit lineaj: relation 
and as an rms scatter of 24% (14% intrinsic). MS1455+22 is plotted as an open triangle in both panels and is excluded in the fits. 



is defined as: 



_ (3- ;)D--r- 



(7) 



where N„et is the background corrected cluster galaxy 
counts, D is the angular diameter distance to the cluster 
redshift, Q is the angular size of the counting aperture, 
is the angular area of the counting aperture, I-y is ge- 
ometric deprojection constant which depends on 7 (I-,, = 
3.78 for 7 ^ 1.8), and is the integrated cluster LF up 
to the apparent magnitude M, which corresponds to an 
absolute magnitude mo at the redshift z. The normaliza- 
tion of the LF (^*) is the universal normalization, and 
not the cluster normalization, so that Bgc is the over- 
density of galaxies compared to the field density, not the 
average cluster density. This formula for B^c is slightly 
different than formula presented in Yee & Lopez-Cruz 
(1999, their equation 3) in that it contains the term. 
We note that in their formula for Bgc, Yee & Lopez-Cruz 
(1999) had inadvertently left out the Aq term due to a 
transcribing error. 

At first, Bgc may appear to be an unnecessarily com- 
plicated measure of the cluster richness; however, using 
it has some distinct advantages. As interest in an obser- 
vationally cheap proxy for cluster mass has increased, so 
have the number of studies which have looked at the cor- 
relation between cluster richness and mass (e.g., Hansen 
et al. 2005; Miller et al. 2005; Gilbank et al. 2004; 
L04; Rines et al. 2004; Ramella et al. 2004; Lin et al. 
2003; Kochanek et al. 2003: YE03). These studies de- 
fine the cluster richness as tlie number of galaxies within 
a fixed aperture, to a fixed magnitude limit. However, 
because different authors use different apertures, and dif- 
ferent magnitude limits, comparison between studies is 
extremely difficult. The advantage of the Bgc parameter 
is that it assumes a universal spatial distribution, and 
luminosity function for clusters. Therefore, aside from 
statistical fluctuations, the value of Bgc is, in princi- 
ple, identical regardless of what aperture is used to count 
galaxies, and which magnitude limit is chosen. 

Yee & Lopez-Cruz (1999) verified this was true, but 
showed that measuring Bgc from a fixed aperture with 
R = 500 kpc (using Hq = 50 km s~^ Mpc~^) produced 
the lowest statistical errors. Furthermore, they showed 
that counting galaxies to M < M^ -|- 1 was all that was 
required for Bgc to be robust. Therefore, our values of 
BgcK are computed using a fixed aperture of 500 kpc 



(although we use Hq = 70 km s~^ Mpc~^) and by count- 
ing galaxies to K* -I- 1. 

Again, because Bgc requires counting the number of 
cluster members, we use statistical background sub- 
traction, rather than the spectroscopic weights. YE03 
demonstrated that Bgc's computed using both techniques 
agree extremely well. In order to determine the inte- 
grated LF parameter (\I>[M(mo,2;)]) a universal cluster 
luminosity function must be assumed. We showed in Pa- 
per I that besides the passive evolution of the stellar pop- 
ulations, this is a reasonable assumption for the CNOCl 
clusters. When calculating ^'[M(mo,-2)] we assume pas- 
sive evolution of the LF of our "average" cluster (K* = 
-24.14, a = -0.84). Unfortimately, we do not have a mea- 
surement of the universal (t>* in imits of Number/Mpc^. 
This cannot be determined from the cluster data alone 
because clusters are high-density regions in the universe, 
and not representative of the mean value of cj)* ■ Rather 
than using values from studies of the K-band field LF 
(which still have fairly large errors) we adopt the ap- 
proach of Yee & Lopez-Cruz (1999) for determining (j>* . 
They showed that a self-consistent value of (j)* could be 
estimated by evolving the cluster LF assuming a sim- 
ple model for the evolution of M*(z). The normaliza- 
tion ((j)*) is determined by requiring that the integrated 
coimts from the LF reproduce the background galaxy 
counts. We perform the same procedure using a z/ =5.0 
model passive evolution model (see Paper I) to parame- 
terize the evolution of M*{z). By comparing the model 
background counts to the true background counts and 
using a x^-maximum-likelihood technique we determine 
(j)* = 4.34 X 10-3 Mpc-^ This value of (j)* is slightly 
larger than the ^* = 3.40 ± 0.29 x 10"^ Mpc'^ deter- 
mined locally (Kochanek et al. 2001), and about 3 times 
as large as the </>* = IJStl'l ^ 10"^ Mpc"^ measured 
at 0.2 < z < 0.65 by Pozzetti et al. (2003). However, 
even if this value is incorrect it will be systematically in- 
correct for all clusters and will only affect the intercept 
in the correlation between Bgc,if and other physical pa- 
rameters. It will not affect the slope or scatter, which 
are of principle interest. 

Column 6 of Table 2 lists the measured values of Bgc 
for the clusters. The errors have been computed using 
the prescription from Yee & Lopez-Cruz (1999), 



ABgcK ^ {Nnet + l-3^jVbg)V^ 

Bgc,K ' 



net 



(8) 
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where Nbg is the number of background counts within 
the 500 kpc and the 1.3^ term is used to approximately 
account for the clustering of background galaxies. 

In the left panel of Figure 5 we plot Log(Bgc,if) vs. 
Log(M2oo) for the clusters. The best fit linear relation 
with MS1455+22 excluded is 

Log{M2oo) = (1.62 ± 0.2i)Log{Bg,,K) + (9.86 ± 0.77). 

(9) 

The fit has a reduced-x^ of 1.29, and the rms scatter 
of the correlation is 35% (18% intrinsic). This implies 
that the K-band selected richness is slightly better than 
the total K-band light for predicting the mass of a galaxy 
cluster. The difference in total and (intrinsic) scatter be- 
tween the two parameters is notable, 35%(18%) in M200 
- Bgc,if vs. 43%(31%) in M200 - ^•200,k and this may be 
because of the smaller aperture used to measure Bgc- 

In the right panel of Figure 5 we plot Log(Bgc,i<:) vs. 
Log(T2,). The best fit linear relation with MS1455+22 
excluded is 

Login) = (0.94±0.13)Lo5(Bgc,K)-(2.20±0.42), (10) 

and has a redueed-x^ of 2.00 and an rms scatter of 25% 
(16% intrinsic). This is very similar to the scatter seen 
in the L2oo,if - T^^ relation and shows that the fixed- 
aperture richness is an excellent indicator of the cluster 
X-ray temperature at the ^ 25% level. 

If we compare the accuracy of BgcK to the optical 
Bgc we find that they are almost identical at predict- 
ing M200 (scatters of 35%, and 31% respectively) as well 
as Tx (scatters of 25% and 21% respectively). Bgc^x 
for the CNOCl clusters also has a similar scatter to the 
fixed- aperture richnesses measured for local clusters in 
the K-band. L04 showed that the scatter in the number 
of galaxies within 0.75 Mpc vs. M500 was 43%, and esti- 
mated that ~ 24% of the scatter was intrinsic. The L04 
cluster masses are measured using and therefore the 
it is more relevant to compare the L04 intrinsic scatter 
to the intrinsic scatter in the CNOCl Bgc - T^; relation. 
The scatter in the L04 relation is somewhat larger than 
the scatter from the CNOCl relation, but this may be 
due to the fact that the CNOCl clusters are generally 
more massive than the L04 clusters. 

These results demonstrate that relatively cheap IR 
imaging can be used to determine cluster masses with 
good accuracy. Specifically, IR imaging of even higher 
redshift clusters {z > 0.5) will be valuable for deter- 
mining masses in large optical cluster surveys. IR rich- 
nesses will likely be more robust than optical richnesses 
at predicting cluster masses at high-redshift because at 
high-redshift optical bandpassess probe bluer parts of a 
galaxy's spectrum where the star-formation properties 
can drastically alter a galaxy's luminosity. 

It is also possible that in the case of large cluster sur- 
veys, Bgc, if may prove to be as valuable for predicting 
the cluster halo mass as the cluster velocity dispersion or 
X-ray temperature. Determining the cluster mass from 
either the line-of-sight velocity dispersion or X-ray tem- 
perature require assumptions about the dynamical state 
of the cluster. Dynamical masses require that the clus- 
ter is in virial equilibrium and that that the shape of 
the velocity ellipsoid is known (generally it is assumed 
to be isotropic). X-ray masses require the assumption of 
hydrostatic equilibrium, and spherical symmetry. They 



also require a correction if the cluster has a cooling flow. 
In cases where these assumptions do not apply (such as 
cluster-cluster mergers, or a collapsing system in the pro- 
cess of forming) the masses of the clusters may be poorly 
estimated, and these clusters will contaminate the mea- 
surement of the cluster mass function. The number of 
catastrophic outliers can have serious consequences on 
the measured cosmological parameters because they are 
particularly sensitive to the number of rare, massive clus- 
ters. On the other hand, if clusters which are unrelaxed, 
or currently undergoing a merger do not have the K-band 
light of their galaxies altered significantly, then using K- 
band richness as a mass indicator could potentially be 
a robust method for determining the masses of clusters 
that are not in dynamical or hydrostatic equilibrium. 

One way to test this hypothesis is to use weak-lensing 
to measure the mass of clusters and compare those 
masses to Bgc, and Bgc.K- Weak-lensing does not re- 
quire any assumptions about the dynamical state of the 
clusters and therefore it is an unbiased (although statis- 
tically noisy) measure of the cluster mass. Weak lensing 
masses for a subsamplc of the CNOCl clusters have been 
measured and will be compared to Bgc^K in a future pa- 
per. 

6. THE K BAND MASS-TO-LIGHT RATIO 

Recent studies of the K-band M/L ratio in local clus- 
ters (Lin et al. 2003, L04, Rines et al. 2004, Ramalla 
et al. 2004) have shown that the K-band cluster M/L 
ratio is an increasing function of cluster mass (although 
Kochanek et al. 2003 find it is roughly constant with 
mass). These studies have also found that the M/L ra- 
tio is a slowly decreasing function of radius, with the 
integrated M/L ratio at R200 (M200/L200.A:) being ^ 
15% smaller than the M/L ratio at R500 (e.g., L04). 
Here we present the first K-band M/L ratios of massive, 
intermediate-redshift clusters. 

We use the k and evolution corrected L2oo,k values to 
compute M2oo/L2oo,A' and therefore it is the M2oo/L2oo,_ff 
ratio of clusters corrected to z — 0. The values are listed 
in Table 3 and the errors computed by propagating the 
M200 and L200 errors in quadrature. In Figure 6 we plot 
Log(M2oo/L2oo.E:) vs Log(M2oo)- There is a clear corre- 
lation of M/L with M200 in the CNOCl clusters. The 
Spearman rank-correlation coefficient for these data is 
0.60, which implies a probability of 0.017 that the vari- 
ables are uncorrelated. The solid line in the plot is the 
best-fit linear relation with MS1455-I-22 excluded: 

Log{M2oo/L2oo,K) = (0.57±0.13)Lo5(M2oo)-(6.92±2.04). 

(11) 

The fit has a reduced-x of 0.964 and implies that M/L 

oc '\/iO-57±o.i3 ^ Interestingly, this slope is about a fac- 
tor of 3 steeper than the M/L oc M°-i^=^° " that would 
be inferred using the Log(M20o) - Log(L2oo,i<:) relation 
(Eqn 5) where the variables are less directly correlated 
(they still both depend on R2oo)- As a comparison, the 
inferred relation is plotted as the dashed line in Figure 
6. The inferred relation, M/L oc M^'^^**^-^^, agrees well 
with the relation from local clusters where M/L cx M" 
with a = 0.26 ± 0.04 (L04), and a = 0.31 ± 0.09 (Rines 
et al. 2004), but does not agree well with the value for 
local groups, a = 0.56 ± 0.05 (Ramella et al. 2004). The 
Ramella et al. relation is similar to the correlation ob- 
tained by fitting Log(M/L) vs. Log(M2oo) directly. 
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3.0 3.2 

Log BgcK 

Fig. 5. — Left Panel: Plot of Log(M2oo) vs. Log(Bgc,K) for the CNOCl clusters. The solid line is the best-fit linear relation and has an 
rms scatter of 35% (18% intrinsic). Right Panel: Plot of Log(Ta;) vs. Log(Bjc,lf ) for the same clusters. The solid line is the best-fit linear 
relation and has an rms scatter of 25% (16% intrinsic). MS1455-t-22 is plotted as an open triangle in both panels and is excluded in the 
fits. 



It is likely that the inconsistent slopes from our own 
data, as well as between the local cluster and group M/L 
vs. M200 arc primarily caused by the large scatter in 
the M200 - 1^200,^ scahng relations and the difficulties 
associated in fitting to data with a larger scatter than 
is accounted for by the error bars. It is probably not 
caused by real differences in the M/L vs. M200 scaling 
relations for the different cluster samples. In our own 
data the M200 - L200,i<r relation has a large rcduccd-x^ 
(2.65). Examination of Figure 4 suggests that the in- 
flated reduced-x^ is not caused by a linear fit being the 
incorrect model for these data, but is caused by a hand- 
ful of significant outliers to the relation. 

Another consideration is that the intrinsic correlation 
between M/L and M200 (mass is in both parameters) can 
contribute to the discrepancy in slope. The Ramella et 
al. (2004) masses and our own masses are determined 
from the cluster velocity dispersion (cri) and equation 2, 
where A4200 a^. The L04 masses are determined from 
using the fit from Finoguenov et al. (2001) where 
M200 oc T^-^^. Given that M200 depends on af but only 
on T^-^^, catastrophic errors in tri will translate to larger 
errors in M200 than catastrophic errors in T^, will. Con- 
sequently, given that the Log(M/L) vs. Log(M2oo) slope 
is < 1, clusters which have their velocity dispersions in- 
correctly measured (and not properly accounted for by 
the errors, such as in the case of non-virialization) will 
steepen the slope of the M/L vs. M200 correlation be- 
cause the mass term is in both parameters. 

As a test, we remove the highest and lowest mass 
CNOCl clusters (MS0440+02 and MS0451-03, both of 
which are significant outliers in the M200 - L200 relation) 
and refit. We find a slope of a = 0.44 ± 0.15, which is 
more consistent with the a = 0.17 ± 0.11 inferred di- 
rectly from the M200 vs. L200.K relation. Unforttmately, 
this problem of outliers makes determining a robust slope 
for the M/L vs. M200 correlation difficult. The data 
are less directly correlated in the case of Log(M2oo) vs. 
Log(L2oo,A:), and therefore we return to that relation and 
adopt M2oo/L200,K o: m"-^^^"-^^ as the best slope for 
M/L vs. M200 scaling relation for the CNOCl clusters. 

This slope, where M/L increases with increasing M200 
supports the scenario proposed by both L04 and Rines et 
al. (2004), where either the star- formation efficiency of 
cluster galaxies is a decreasing function of cluster mass, 
or else that the amount of intra-cluster light is an increas- 
ing function of cluster mass. Interestingly, the slope of 
this scaling relation for the CNOCl clusters is consistent 



with the the slope in local clusters and this shows that 
the M/L vs. M200 correlation is in place by at least z 
0.3, roughly 4 Gyr in lookback time. 

7. Qm from the oort technique 

The original purpose of the CNOCl project was to 
measure the cosmological density parameter flrn,a = 
Pol Pc using the Oort (1958) technique. The cluster M/L 
ratio, divided by the M/L ratio for closure (M/L)c = 
Pc/i, where j is the field luminosity density, provides 
a direct meastire of Jlm.o which is independent of the 
Hubble parameter. There has been some concern with 
this method because recent numerical simulations sug- 
gest that light is a biased tracer of dark matter, and that 
the M/L ratio of any region in the universe which con- 
tains galaxies automatically provides an underestimate of 
the universal M/L ratio, and an underestimate of ^m,o 
(e.g., Ostriker et al. 2003). However, this is only a sig- 
nificant concern for low-density regions in the universe 
(i.e., the group scale or less). Unlike those lower den- 
sity regions, rich clusters arc assembled from regions of 
> 10 Mpc across (Carlbcrg et al. 1997) and therefore, 
provided that observations extend to the entire volume 
of the cluster (i.e., R ~ R200)) clusters should contain 
a sufficient collapsed volume to provide a representative 
sample of the universal M/L ratio (e.g. Carlbcrg et al. 
1996, Carlbcrg et al. 1997). This is yet to be verified 
in numerical simulations as, until recently, they lacked 
enough volume to contain rich clusters. The CNOCl 
project was designed specifically with this biasing con- 
cern in mind and therefore is a wide-field study of a sam- 
ple of rich clusters. 

Although biasing is not a major concern for the clus- 
ter technique, we know that the stellar populations of 
galaxies in the cluster and field environments are dif- 
ferent (e.g., EUingson et al. 2001, Balogh et al. 1999, 
Poggianti et al. 1999, Dressier et al. 1997) and without 
a good understanding on how these differences affect the 
integrated luminosity of these regions, there arises the 
potential for a systematic error in 0.mfi- 

Using the r-band data for the CNOCl clusters, Carl- 
bcrg et al. (1997) showed that n,„.o = 0.19 ± 0.06 (ran- 
dom) ±0.04 (systematic). In their analysis they found 
that the r-band cluster galaxies were on average 0.11 ± 
0.05 mag fainter than field galaxies at this redshift and 
hence they decreased the cluster M/L ratio accordingly. 
Given the need for this correction, it is advantageous to 
perform the same analysis with infrared data, because 
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Fig. 6. — Plot of Log(M2oo/L200,if ) vs Log(M2oo) for the CNOCl clusters. The solid line is the best-fit relation. The dashed line is the 
relation that is inferred from the Log(M2oo) - Log(L200,ir) fit. MS1455+22 is plotted as an open triangle and is excluded in the fit. 



infrared light is a better tracer of total stellar mass than 

optical light. The infrared luminosity of galaxies docs 
not depend strongly on the current star-formation prop- 
erties of the stellar population and consequently, little or 
no correction is required to account for the varying star- 
formation properties of the field and cluster environment. 
The only potential correction required would be if the 
cluster environment played a role in the star-formation 
efficiency of galaxies. A differing star-formation effi- 
ciency could manifest itself in two ways: 1) the stellar 
mass may be distributed differently in the cluster /field 
environment (i.e., the shape of the field and cluster LFs 
may be different); or 2) baryons may be converted into 
stars at a different rate than in the field (i.e., the ratio of 
the stellar mass to dark matter mass could be different 
between these environments). 

In terms of problem 1), Paper I showed that the cluster 
K* was brighter than, but consistent with the field value 
(AK* = 0.25 ± 0.26 mag). Paper I also showed that a 
from the CNOCl cluster LF was the same as a for local 
clusters as well as the local field. Those results agreed 
well with local measurements which showed the cluster 
K* and a are roughly the same as the field (although we 
note that both L04 and Rines et al. 2004 find that the 
cluster K* is ~ 0.2 mag brighter at the 2(7 level). Given 
that both K* and a for the cluster LF are similar to the 
field LF, the distribution of galaxies in terms of luminos- 
ity (and by corollary stellar mass) should be similar in 
both environments. 

In terms of problem 2), it is important to note that 
the cluster M/L ratio continues to be an increasing func- 
tion of cluster mass, even for objects as massive as the 
CNOCl clusters. This shows that the overall rate of 
converting baryons to stars is a function of environment. 
The rich CNOCl clusters have a K-band M/L ratio which 
is about an order of magnitude larger than the M/L ra- 
tio of groups (e.g., L04, Ramella et al. 2004), and in 



turn, the group M/L ratio is about an order of magni- 
tude larger than for individual galaxies (e.g., Brinchmann 
& Ellis 2000). Remarkably, despite the fact that rich 
clusters are two orders of magnitude less efficient in con- 
verting baryons to stars than individual galaxies, their 
continually increasing M/L ratio with increasing mass 
suggests the possibility that they may still be slightly bi- 
ased tracers of stellar mass in the universe. Unless the 
M/L vs. M200 relation flattens at the scale of ^ 10^'"' 
M0, the M/L ratio of clusters could be lower than the 
universal value. Unfortunately, there is no data for the 
K-band M/L ratio of superclusters to determine if the 
M/L ratio flattens at this size scale. 

The mean M/L ratio for the clusters (corrected for pas- 
sive evolution to z = 0) is 61.2 ± 1.8 M0/L0. Combin- 
ing this with the luminosity density of the local universe 
measured by Kochanek et al. (2001) we find 0^,0 — 
0.22 ± 0.02 (random). This result is in good agreement 
with combined constraints from WMAP third-year re- 
sults (WMAP3) and various other techniques (Spergel 
et al. 2006). The combined WMAP3 constraints on 
fim.a range from a low-density of Qm.a 



0.226tro« to 



a high-density of flm.o = 0.299_q q25. Interestingly, the 
measurement from the clusters is closest to the lowest 
density of the of preferred values from combined con- 
straints. Previous studies using K-band M/L ratios of 
local clusters have also found values of 0^,0 on the low- 
end of preferred values. 

Both Rines et al. (2004) and Lin et al. (2003) mea- 
sured Qm,o using their local cluster samples. Rines et al. 
found that the result depended strongly on the location 
in the clusters where the M/L ratio was measured. The 
cluster infall region had a lower M/L ratio than the vi- 
ralized region and therefore they found flm,o = 0.18 ± 
0.04 from the virial region and flm,o = 0.13 ± 0.03 from 
the infall region. Lin et al. (2003) performed the same 
analysis with a larger sample of 27 clusters. They found 
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f^m.o = 0.17 ± 0.02 using the mean M/L ratio of all clus- 
ters, and rimfl = 0.19 ± 0.03 using a subsample of their 
most massive clusters. 

Our analysis is in good agreement with both these 
studies, as well as the r-band measurement from the 
same clusters by Carlberg et al. (1997). Our results 
further confirm that Qm.o from clusters agrees well with 
concordance values but tends to be on the low-density 
end of preferred values. It is possible that the cluster 
measurements prefer lower ri^.o because their M/L ra- 
tios are still lower than the universal value. A measure- 
ment of the K-band M/L ratio on the supcrcluster scale 
would be a useful way to test this hypothesis. Regard- 
less, the fact that the cluster ilm,o still agrees well with 
recent measurements of ilm,o using a variety of indepen- 
dent techniques, and combinations thereof suggests that 
the cluster M/L ratio is unlikely to be significantly lower 
than the universal value. 

8. SUMMARY 

We have presented the K-band scaling relations for 15 

modcratc-redshift clusters with extensive optical spec- 
troscopy and wide-field K-band imaging. The cluster 
HON is well-correlated with M500 and T^,; however, the 
intrinsic scatter in the scaling relations at z ^ 0.3 is fairly 
large (37% and 46% respectively). Comparing to local 
clusters we find that the HON is consistent with no evo- 
lution at fixed cluster mass out to z ~ 0.5. This result, 
in tandem with the purely passive evolution of the clus- 
ter LF, and the fact that the cluster LF does not depend 
on mass suggests that if the significant tidal stripping 
of galaxy halos in clusters seen in numerical simulation 
occurs, the stellar mass within the halo is tightly bound 
and remains intact. This interpretation is further sup- 
ported by recent SPH simulations which show that the 
baryonic matter is the amongst the most tightly bound 
mass within the halo. 

Our data also show that both L2oo.if and Bgc.K arc 
well-correlated with the cluster dynamical mass and X- 
ray temperature. The slope of the L200 - M200 relation at 
z ~ 0.3 (Ij2oo,k oc M2oq'^^ '^^) is consistent with the slope 
measured for local clusters, suggesting that the cluster 
scaling relations are in place by at least z ^ 0.3. The 
good correlation, and relatively small scatter (intrinsic 
-|- measurement) in the Bgc^x - M200 relation (35%) and 
Bgc,K - T^x relation (25%) suggests that Bgc^K would 
be useful mass indicator for upcoming cluster cosmol- 
ogy projects. 

Our results from the M/L ratios of the clusters show 
that the moderate-redshift CNOCl clusters are very sim- 
ilar to local clusters in that the M/L ratio is a slowly 
increasing function of cluster mass. By comparing the 
cluster M/L ratio to the local luminosity density we es- 
timate that flm.Q = 0.22 ± 0.02, which agrees well with 
the original analysis of the CNOCl clusters using optical 
data as well as with the local estimates using K-band 
photometry The measured value is also in good agree- 
ment with recent WMAP third year joint constraints; 
however, similar to previous cluster studies the value is 
on the low-density end of preferred values. 

The combined analysis of this paper and Paper I 
present a relatively simple picture of the evolution of the 
near-infrared properties of clusters from z = to z ~ 0.3. 
The correlation between the cluster near-infrared prop- 



erties (i.e., L2oo,K, N500, and M2oo/L2oo,k) and M200 
shows no significant change between z = and z ^ 0.3. 
Furthermore, the scatter in these scaling relations is sim- 
ilar at both redshifts. The cluster LF shows only passive 
evolution between z = and z ^ 0.3 and does not de- 
pend on cluster mass. In addition to this, the (small) 
difference between the field and cluster LF at z = is 
unchanged at z ^ 0.3. These results all show that there 
is little evolution in the the bulk of the stellar mass in 
cluster galaxies over this redshift range besides the pas- 
sive aging of the stellar populations. 

Specifically, it appears that 1) the significant tidal 
stripping of halo in high-density regions seen in N-body 
simulations does not affect the stellar mass contained 
in galaxies nor change the cluster scaling relations; 2) 
given passive evolution of the LF, and no-evolution in 
the HON, mergers and disruptions arc unlikely to play 
a significant role in cluster galaxy evolution at z < 0.3; 
and 3) that the changes seen in the cluster stellar popula- 
tions over the redshift interval z = to z ^ 0.3 (i.e., the 
increase in blue- fraction/star- formation properties, and 
evolution of the morphology-density relation) are "su- 
perficial" - they result from changes in a small part of 
a galaxy's stellar mass, while the majority of the stellar 
mass is already in place and passively evolving. 

Overall, it appears that the bulk of the stellar mass in 
cluster galaxies is in place, and evolving passively with 
few mergers or disruptions between z ^ 0.3 and z = 0, 
and that the the cluster scaling relations are produced 
by processes that occur at higher redshifts. 
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TABLE 1 

Physical Properties of the CNOCl Clusters 



Cluster 


z 


km 


R200 
Mpc 


M200 
(5) 


keV 


A2390 


0.2279 


1095 ± 61 


2.00 


20.1 ± 


3.4 


10 3+"-® 


MS0016+16 


0..5466 


1243 ± 128 


1.89 


24.0 ± 


7.5 


9 2+1° 


MS0302+16 


0.4246 


656 ± 93 


1.07 


3.7 ± 


1.6 


4+3 


MS0440+02 


0.1965 


611 ± 62 


1.13 


3.5 ± 


1.1 




MS0451+02 


0.2010 


979 ± 76 


1.81 


15.0 ± 


3.5 




MS0451-03 


0.5392 


1354 ± 105 


2.06 


33.3 ± 


7.6 


10.211:° 


MS0839+29 


0.1928 


788 ± 104 


1.46 


6.6 ± 


2.7 


4+0.3 
^-0.3 


MS1006+12 


0.2605 


912 ± 101 


1.63 


11.2 ± 


3.8 




MS1008-12 


0.3062 


1059 ± 107 


1.85 


17.2 ± 


5.3 




MS1224+20 


0.3255 


798 ± 90 


1.38 


7.5 ± 


2.5 


A q+0.6S 


MS1231+15 


0.2350 


662 ± 69 


1.20 


4.5 ± 


1.4 




MS1358+62 


0.3290 


910 ± 54 


1.57 


12.0 ± 


2.1 




MS1455+22 


0.2570 


1169 ± 140 


2.10 


22.0 ± 


8.2 


4.4±- 


MS1512+36 


0.3726 


697 ± 96 


1.17 


4.7 ± 


2.0 


3.41°:? 


MS1621+26 


0.4274 


833 it 55 


1.36 


6.8 ± 


1.4 


7+3 
'-2 
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TABLE 2 

NIR Properties of the CNOCl Clusters 



Cluster 


T 


^ T 

e 1^200, K 


(M/L)200,if 




Hgc,K 




N500 


e N500 




T -1 n1 ^ 
LqX 10" 


LqX lO^'^ 






Mpc^° 


Mpc^-° 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 




3.14 


0.25 


64 


14 


2302 


303 


338 


50 




3.94 


0.43 


61 


24 


2298 


311 


389 


56 




0.81 


0.17 


45 


25 


674 


157 


63 


25 


Msn44n-i-n2 


1.63 


0.43 


21 


8 


1083 


210 


125 


25 


MS0451+02 


2.40 


0.29 


63 


18 


1571 


252 


313 


44 


MS0451-03 


2.79 


0.08 


119 


34 


2897 


343 


305 


59 


MS0839+29 


0.79 


0.14 


83 


43 


1768 


266 


207 


34 


MS1006+12 


1.96 


0.18 


57 


24 


2496 


316 


377 


43 


MS1008-12 


2.64 


0.28 


65 


28 


1745 


266 


198 


44 


MS1224+20 


1.34 


0.27 


56 


24 


733 


169 


139 


32 


MS1231+15 


1.23 


0.14 


37 


14 


967 


198 


130 


27 


MS1358+62 


3.05 


0.25 


39 


8 


1866 


275 


269 


41 


MS1455+22 


2.03 


0.37 


108 


50 


982 


199 


89 


46 


MS1512+36 


0.65 


0.10 


71 


37 


772 


172 


49 


25 


MS1621^26 


2.17 


0.15 


27 


7 


154.'^ 


219 


225 


38 
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